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We consider a system of two solid state charge qubits, coupled to a single read-out device, con- 
sisting of a single-electron transistor (SET). The conductance of each tunnel junction is influenced 
by its neighboring qubit, and thus the current through the transistor is determined by the qubits' 
state. The full counting statistics of the electrons passing the transistor is calculated, and we dis- 
cuss qubit dephasing, as well as the quantum efficiency of the readout. The current measurement is 
then compared to readout using real-time detection of the SET island's charge state. For the latter 
method we show that the quantum efficiency is always unity. Comparing the two methods a simple 
geometrical interpretation of the quantum efficiency of the current measurement appears. Finally, 
we note that full quantum efficiency in some cases can be achieved measuring the average charge of 
the SET island, in addition to the average current. 
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I. INTRODUCTION 

A solid state charge qubit, formed by a single electron 
trapped in a double quantum dot (DQD), is an interest- 
ing system for studying quantum coherence in the solid 
stateii^i^ii. 

The state of charge qubits can be read out by sensitive 
electrometers, such as the radio- frequency single-electron 
transistor (RF-SET)-. The straightforward scheme is to 
couple the charge qubit capacitively to the small SET 
island^. The SET acts as a variable resistor, dependent 
on the qubit state. The SET is included in an LC-tank 
circuit, and the qubit state can be detected by measur- 
ing the tank circuits dissipation. Theoretical estimates 
for the back-action on the qubit indicate that single-shot 
qubit readout is possible^'®. 

In principle the RF-SET can also be used to detect 
the charge of a QD^"^''. But a conceptually even simpler 
method is to place the QD close to a quantum point 
contact (QPC)ii. The charge of the QD will influence 
the transparency of the QPC, and by applying a driving 
voltage across the QPC, the charge can be determined 
through a current measurement. This method has been 
successfully used to read out the charge of a DQD^ii^i^. 

Starting with Ref. [ij, there are numerous theoretical 
works on the measurement induced dephasing rate of 
a charge qubit read out by a QPC, see e.g. Ref. [H and 
references therein. The dephasing rate can be compared 
with the measurement time tms, and from fundamental 
principles of quantum measuremeniji^, the quantum ef- 
ficiency 7] = {tms^ip)^^ has an upper bound of unity, 
implying that one cannot distinguish two states without 
destroying the quantum coherence between them. 

An interesting combination of the SET and QPC mea- 
surement techniques was proposed by Tanamoto and 
Hui^i, and is also the focus of our present paper. The 
setup allows for reading out two charge qubits using a 
single SET. The two qubits are positioned close to the 
SET source and drain junction respectively (see Fig.[T]). 
Thus the state of the left/right qubit will infiuence the 
conductance of the left/right timnel junction. By ap- 



plying a source-drain voltage and measuring the current, 
both qubits can be read out. In Ref. 17 the focus was 
made on the transient dynamics of qubit-detector system, 
as well as the ensemble averaged current. 

In this paper we concentrate on the single-shot read- 
out properties, by calculating the measurement times 
from the noise and full counting statistics^* of the charge 
transport. The quantum efficiency is obtained by com- 
paring with the measurement induced dephasing rates. 
From the quantum inefficiency, in detecting certain qubit 
states, we conclude that the qubits get entangled with de- 
grees of freedom not measured in current measurements. 
This leads to an analysis of a measurement setup, where 
instead of the current, the charge of the SET island is 
detected, in real-time. We find that this measurement 
always give full quantum efficiency. 

The outline of the paper is as follows. In Sees. HIl 
and mil we derive a master equation for the dynamics of 
the density matrix of the detector-qubit system, using 
an open quantum systems approach^^. In Sec. llVl we use 
this formulation to obtain the single-shot measurement 
time, as well as the dephasing rates, and consequently 
also the quantum efficiency. We then apply the general 
results to a few analytically tractable cases. In SecfVlthe 
analysis is extended to the case when the charge of the 
SET island is monitored. The properties of that setup 
is compared with the current measurements in Sec. IVIl 
leaving the conclusion and discussion for Sec. IVIII 



II. THE MODEL 

We consider two double-dot charge qubits, read out by 
a quantum dot single-electron transistor. The SET is- 
land contains one spin-degenerate level with energy E^, 
and the Coulomb interaction gives rise to an additional 
energy U when it is doubly occupied. Each qubit is elec- 
trostatically coupled to one of the SET's quantum point 
contacts, as shown in Figure [T] There is also an elec- 
trostatic interaction between the qubits, causing a direct 
qubit-qubit coupling with strength J. The two qubits 
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FIG. 1: A Coulomb blockade island in a single-electron tran- 
sistor geometry. The conductance of each tunnel junction is 
influenced by a nearby charge qubit. (See text.) 

have a Hamiltonian of the following form 

Hqb— ^ {^a^yax + ^a<^az) + Jl^Lz(^Rz, (1) 
a=L,R 

where Uax and CTq-z are Pauli matrices (a € f^Q. 
denotes the tunnel coupling strength between the two 
dots of each qubit, and Aq, the charging energy difference 
within each qubit. We can also express the Pauli matrices 
in terms of the creation and destruction operators of the 
extra electron in the qubit quantum dots, 

CTax = a\ba + b\aa, CTaz = " ^a^a, (2) 

where a\ and aa denote the creation and destruction 
operators in the upper dot of the a S {L, R\ qubit, and 
6jj and ha operates on the lower dot. 

The SET Hamiltonian can be split into three pieces: 
HsET — Hres + His + Ht ■ Hres dcscribcs the left and 
right reservoirs 

a=L,R;s=T,l k 

where Eka is the electron energy in the a reservoir at 
wave vector k, and cj.^^ and Ckas are corresponding elec- 
tron creation and destruction operators. For the SET 
island, we have, 

H,s = ^ddlds + Ud\d^d\di, (4) 

s=^,^ 



where d\ and ds are creation and destruction operators 
of the electrons on the island. We consider the transport 
through the QPCs to be in the low transparency regime, 
described by a tunnel Hamiltonian, 

a=L.fi;s=T,i fe 

(5) 

where we assume that the tunnelling strengths Vka of 
electrons between the reservoirs and the Coulomb island 
are independent of the spin degree of freedom, and the 
spin s is conserved during the tunnelling processes. The 
counting fields Xa for a G {L, R} are introduced to keep 
track of the number of electrons mj^j^ which passed 
through the LjR junction from the SET island to the 
LjR lead. The operators e^*-^" changes this quantum 
number by one, e^^^" IrUa) — \ma ± 1). 

The qubits electrostatically interact with the QPCs, 
and influence the tunnelling rates according to the fol- 
lowing Hamiltonian, 

a,s,k 

where SVka are coupling strengths between the QPCs 
and the qubits. Depending on the qubits' states, the 
tunnelling strengths of the left and right QPCs change 
from Vka to Vka ± SVka ■ For simplicity, we neglect the 
relative phase between these tunnelling strengths, but we 
still allow for a weak energy dependence of their magni- 
tudes. 



III. A KINETIC EQUATION 

One conventional approach to continuous quan- 
tum measurement problems is called the open system 
approach^. In general terms, an open system is a lim- 
ited quantum system coupled to another quantum sys- 
tem with a large number of degrees of freedom, called 
the environment. In our case, the reservoirs (leads) con- 
nected to the left and right QPCs constitute an environ- 
ment for the limited quantum system consisting of the 
two qubits and the SET island. The system hamiltonian 
is Hs = Hqb + His , while the hamiltonian of the environ- 
ment is Hres and the interaction hamiltonian is given by 
the tunnelling terms Hint — Ht + H^ . 
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The starting point is the Liouville equation on integro-difFerential form 



dtPtotit)=-i[H,„tit),ptot{0)]- / dt' H,r,t{t),[H,ntit'),ptot{t')] 

Jo L 



(7) 



where Hint{t) is the interaction hamihonian and ptot{t) the density matrix in the interaction picture. Then we 
proceed to a markofhan kinetic equation for the reduced density matrix of the system p{t), making the Born and 
Markov approximations, 



dt'Tu, 



Hint{t), [Hi, 



tin Pit)]]}. 



(8) 



where Tr^es indicates taking the trace over the reservoir degrees of freedom, which are taken to be in thermal 
equihbrium, at their respective chemical potential. The integrals on the right-hand side of Eq. ([5]) can be approximately 
evaluated, taking the lower limit of the integral to minus infinity. They have in general both real and imaginary parts, 
where the imaginary parts effectively can be absorbed into renormalized system energies. The real parts give the 
dissipative dynamics, i.e. the electron tunnelling rates. A typical example of a term contributing to the tunnelling 
rate from lead a to the island is 



(9) 



where we have taken the sum to an integral by introducing the density of states in the leads ga (E) , and the energy 
dependent tunnelHng amplitude Va{E), and fa{E,T) is the Fermi function of lead a. The bias voltage between the 
left and right lead enters through the shifted Fermi distributions. For details, please see Ref. [20l . 



A. The non-zero tunnelling rates 



Since our quantum system has 16 degrees of freedom the reduced density matrix has 256 elements, and thus the 
kinetic equation in general 256 x 256 terms, we have to make some simplifying assumptions to arrive at a tractable 
set of equations. 

First we consider a source-drain voltage across the SET large enough compared to temperature so that we may 
neglect backwards tunnelling. Thus, without the coupling to the qubits, and with spin-independent tunnelling, there 
are four different rates to consider; tunnelling onto the unoccupied SET island from the left lead (Fl), tunnelling onto 
the singly occupied SET island from the left lead (F^), tunnelling off the singly occupied SET island to the right lead 
{Tjij, and finally tunnelling off the doubly occupied SET island to the right lead (F'^^), see Fig. [5] 
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FIG. 2: The four non-zero tunneling rates for electrons, on and off the SET island. The applied source-drain voltage eV creates 
a difference in chemical potentials between the left and right lead eV = fiL — pr- 

Including the coupling to the qubits each rate will acquire a state-dependent shift Fj — Fq±5Fq and F^''= — F^±(5F'^, 
see Fig. [21 Not to complicate the system further we assume that qubit energies Aa,i^a, and J are small enough so 
that the tunnelling rates including a flip of the qubits's state are identical to the ones without flipping the qubit state. 

We are now ready to present the kinetic equation for the elements of the reduced density matrix. As basis we chose 
the four qubit product states \A) = \ ||), \B) = \ |t), \C) = \ ti), \D) = \ tt), the SET island states corresponding to 



FIG. 3: An illustration of the capacitative coupling between a double quantum dot qubit and low transparency point contact. 
Left: When the electron is located in the dot closest to the PC, the barrier is higher, leading to a lower tunnelling rate 
r~ = r — sr. Right: For the qubit in the other state, the barrier is lower and the tunnelling rate higher F"*" = F + ST. 



no electron |a), one spin- up/down electron and doubly occupied |c), and the number of electrons which has 

tunnelled across the right junction \m): 



P"iz2(™) = (™l ® (^1 ® (-^il P 1-^2) ® \u) ® |m), 



(10) 



where 21,2:2 G {A, B,C, D} and u G {a, 6|,6|,c}. Here we only keep terms off-diagonal in the qubit states. It is 
easy to verify that terms off-diagonal in the island state do not couple to the diagonal terms, and they decay on a 
time-scale given by T^^^. Because the kinetic equation is translationally invariant in m-space, we Fourier transform 

it with respect to this variable ^"^^^(fc) = J2m e'^™p"^^^ (to)^ '^"'^^^^ , and arrive at the following kinetic equation for the 
elements of the reduced density matrix 



where 



P2122 



^ L 



~^^L [Pgi{zi),Z2 " P''zi,g,{z2) 

^{'^Z2 Jzi ) 

i{Jz2 



ygA.z\),Z2 I^zi,gr(z2) 



zi-pZ2 ik ( 



Z2 + PziZ2 



Pz\z2 - i^R 



Pgr(zi),Z2 Pzi,gr{z2) 



Pgiizi), 



Jz 



Pz 



^ L 



Wzi-p'z2 ik c 

^ R ^ R PziZ2^ 



^Z2 



T^Z2 

^ R 



PziZ2 



bl bl 

Pgr{zi),Z2 Pzi,gr(z2) 



Pgi(zi)-Z2 Pzi,gi(z2) 



n^'LPl,Z2 + ^^R'^R'e^'pl,Z2^ 



V{Jz2-Jz,)~{T^'+V^^)]pl,.^ 



Pgr(zi),Z2 Pzi,gr(z2) 



Pgi{zi),Z2 Pzi,gi{z2) 



^■Z,^Z2 



{pl\ 



Z2 + PziZ2 



(11) 



(12) 



(13) 



(14) 



= r' 



L ' 



^ R. 



r+ 



Ja = -Al - Afl + J, Jb = -Al + Afl - J, Jc - Al - A^, -J, Jd ^ Al + Ar + J, (15) 

and 

gi{A)^C, gi{B) = D, giiC) = A, gi{D) = B, gM)^B, gAB) = A, gAC) = D, gAD) = C. (16) 

We note that, neglecting the counting field (fc ^ 0), the above equations coincide with the rate equations derived by 
Tanamoto and Hu (see Appendix in Ref. ,17,) , using a different approach^"'. Eqs. pip to form one of the main 
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results of this paper, and they constitute the basis of our further analysis. One may note that although the rates 
r^j,^ and ^r^/^ are intimately connected to the microscopic tunnelling hamiltonians in Eqs. ([5]) and ([6]), we may 
now consider different operating regimes in terms of the value of these rates themselves, rather than specifying the 
coefficients of the tunnel hamiltonians. 



B. Numerical methods 



To solve for the time-dependence of the fc-dependent reduced density matrix, we rewrite the kinetic equations on 
matrix form p = M ■ p as follows 



( p\a \ 

Pab 

\pddJ 



-2T-, 



iflR 2i{AH - J) - 2TI 



64 



f pU\ 

Pab 
\pdd) 



V ... -2r + /,,^,, 

Thus the reduced density matrix elements at arbitrary time t are given by the evolution matrix M 



(17) 



(18) 



By considering the time-derivative of the average number of electrons which have tunnelled across the right junction 
(to) — Tr{'m\'m){m\p} = Tr{—idkp{k = 0)} we get the expression for the time-dependent average current 



(Init)) = e^{n 



pim+pim 



(19) 



which also agrees with the one used in Ref. Il7|. To discuss the single-shot measurement time, needed to separate differ- 
ent qubit states, we will also need the time-dependent probability distribution of the number of electrons transferred 
through the detector P{m, t). This we obtain by tracing out the qubits' and island's degrees of freedom. 



P(m, t) = p{z, z; M, u; to, to; t), 



(20) 



which is straightforward using the Fourier transformed fc-dependent density matrix, 

dk 



P{m,t) 



(21) 



We are now ready to discuss the properties of reading out two charge qubits using a single-electron transistor, in a 
few different parameter regimes. 



IV. READ-OUT IN THE QUBITS' EIGENBASIS 

To approach a quantum non-demolition (QND) measurement for a weakly coupled read-out device, it is necessary 
to measure in the qubits' eigenbasis. For the setup considered here, the measurement is performed in the z-direction 
of each qubit, i.e. the measurement determines in which dot the electron of each qubit is located. Thus we consider 
the case when the inter-dot coupling is switched off, i.e. = = 0. We then start from Eqs. (|1HI16P and derive 
simplified kinetic equations. Due to the spin degeneracy of the tunneling, we may treat the two singly occupied island 
states together by defining p\^z2 ~ ^^1^2 + P^Az2- Also, we note that the terms proportional to «(Jz2 ~ Jz-i) vanish by 
going to the rotating frame, with respect to these coherent z-rotations, i.e. PziZ2{t) = e~*^''*2-'^2i)*p^^^^(t). We arrive 
at the following equations for the density matrix in the rotating frame p(t), 



= - 2 + '^^n^LPl... + 2Vr«"r -p^^,^, (23) 



pl,.. = -(r^^ + ^'r')pI^.. + ^^'l'^'l'p\...^ (24) 



where we have suppressed the 'tilde' for notational brevity. 
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A. Dephasing rates and the quasistationary distribution 



Since the interdot coupling is switched off (Aj^ = = 0), the equations for density matrix elements with different 
qubit state indices (212:2) separate. To analyze these equations further we now consider the elements of the reduced 
qubit density matrix PziZ2 — PtiZ2 + ^2122 pIxz^^ obeying the equation 



r ^1^2 'in'' 
(pL JrziZ2 



ipR I Z1Z2 ' 



where we defined the basic dephasing rates 



-pZi 

^ L/R 



n^/B) and r = 



y' Z2 
^ L/R 



(25) 



(26) 



It is clear that these rates are only non-zero for off-diagonal elements (21 7^ 22). Thus the diagonal elements pzz are 
conserved, i.e. p'zz = 0, as required by probability conservation. However, all off-diagonal elements decay, which is 
the definition of dephasing. From this we may conclude that the measurement will eventually destroy all coherence 
between the different qubit product states | J,|), | |t)i I Ti) and | tt)- The expressions in Eq. (pS]) are similar to the 
ones presented in Ref. llTI . and for each junction and tunnel event separately the rates agrees with the result for a 
single qubit read out by a single QPC>^ 

In the regime of weak measurement j^/^/T j^/^ ^ I, the dephasing rates are small 



r 



i + o 



r2 



(27) 



so the dynamics will have two different timescales. On the short timescale given by r^^jj, all reduced density matrix 
elements will be approximately conserved, and we obtain the quasistationary distribution 



p'i,z2{t) 




^ R ^ R 




PziZ2{t) 


-pZlZ2Y' Z1Z2 


|_2p2l22p^l22 _ 


R R 


Pl,Z2it) 




0"p2i2;2-p 2122 

L ^ R 




Pz^Z2 (*) 


T^Z\Z2J-.' Z\Z2 

^ L ^ L 


1- 2r^i^^r^i^^ - 


I -pZiZ2W Z1Z2 

^ ^ R. ^ R. 






p2i22p'2i22 

^ L ^ L 




PziZ2it) 


p2l22p'2i22 


|_ 2p2l22p^l22 _ 


1 T^ZiZ2W Z1Z2 

^ ^ R ^ R 



(28) 



where r--=r-^ 



_-p22 p 2122 "p 21 

L/R^ ^ L/R ~ ^ L/R' 



P 22 

L/i?,' 



, and P"j, is the probability to find the island in the state a S {a, 5, c} 

given that the qubits are in the product state [2) e {|j4), |_B), |C), |Z))}. In order to get a qualitative understanding it 
is useful to note that According to Eq. ((25| . the dephasing will take place on the much longer timescale r^/j^/Jr^^^, 
and as we will see, this is also the time-scale for the measurement. Inserting this quasi-stationary distribution into 
Eq. ([^5]) we find the total dephasing rate between the product states \zi) and [22), 



pZl22 _ _ 
^ 2 



pa 2r^i5' 

Z1Z2 '^L 



Pi 



-pZiZ2 _|_ "p Z\Z2 \ 1 pC O-p Z1Z2 

^ ipR ^ ^ tpL I ~^ ^ziZ2^^ ipR. 



(29) 



One simple interpretation of this formula is that the coherence between state 21 and 22 will decrease only due to 
tunnel events across a junction where the two states differ. 

B. Average current, noise and the measurement time 

The average current will approach a steady-state value on the short time-scale given by T"^^ (see Fig. 2]). In the 
case when the qubits initially are in the product state \z) the steady-state average current 



Pzz^R 



PL2V 



R 



(30) 



is obtained by setting pzz (t) = 1 in Eq. (pS)) and then inserting p"^ , p' 



b 

zz ' t^zz 



p\i and p% in Eq. (HH). 
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FIG. 4: Left: The average current through the SET for the different initial qubit product states, with no qubit interdot coupling 
(A = 0), and SFl ~ SFr = 0.2r. Right: The corresponding probability distributions P{m,t) for the number of electrons m 
having passed the SET at time t — 60/r. 

For times much longer than T^^^, when the average number of electrons which has tunnelled through the SET is 
large (m) ^ 1, the distribution P{m, t) will approach a Gaussian with average (m) ~ ht/e, and a width given by the 
zero-frequency current noise Si. The current noise is readily obtained from the full counting statistics as^^^ 

Si = 2e2 (^_^Ao(fc)) , (31) 

where Ao(fc) is the eigenvalue of the evolution matrix M which approaches zero in the limit k 0. The noise is 
sometimes expressed in terms of the Fano factor /, through the relation Si — 2elf. After some algebra we obtain an 
analytic expression for the Fano factor corresponding to the qubit product state \z), 

sf i^R + r/)' (in? + i^'i?) - (nn + + 4(ri)^(rj)2 

fz = — = — ^ . (32) 

We now define the time needed to distinguish between two different product states \zi) and \z2) as^i 

after which the two corresponding counting distributions are distinguishable. Distinguishable implies that the overlap 
of the counting distributions fulfills^''' 

^ v/^zi(m,i^f )P,,(m,aif ) < -, (34) 



e 



where Pz{m,t) corresponds to the qubits' state \z). 

According to fundamental principles of quantum mechanics, two different states cannot be distinguished classi- 
cally before the quantum coherence between these states has decayedJ^ With our definitions this inequality reads 
r^^^^^mJ^ > li leading to the definition of the quantum efficiency 

V^'^' = pii^ < 1' (35) 

for distinguishing the two states \zi) and \z2), using our suggested measurement scheme. One may note that other 
definitions of distinguishability, and thus measurement time, gives other numerical values for the limit of the quantum 
efficiency. This is the reason why the quantum limit sometimes in literature is given as 1/2. The foundation is now 
established for discussing the properties of the suggested readout scheme in some analytically tractable, and relevant 
cases. I 

C. A single junction detecting a single qubit obtained from our system if the right qubit is uncoupled 



For comparison wc may check the limit of a single QPC 
measuring the state of a single qubit. This situation is 
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{STji = 0), and if the escape rate through the right junc- 
tion is much larger than the in-rates F/j ^ F^. To lead- 
ing (zeroth) order in Vf^/Tji the island is always empty, 
i.e. = 1, the average current is 1^ = 2eF|^, and 

the Fano factor is one. The time needed to separate the 
product state which differ in the state of the left qubit is 
obtained from Eq. (l33l) . 



AC 



t. 



AD 



BC 



t 



BD 



(36) 



5Tl 



l + O 



where the last equality is relevant in the weak coupling 
regime (^F^ <C F^). The corresponding dcphasing rates 
are given by Eq. (pS)) . 

AC 



(37) 



2 Fr - 



6Tl 



l + O 



Thus we find that this readout is indeed quantum limited 
: ri^ c — rf^^ = Tj^^ = 1^ which was first discussed 



AC 



in Ref. [IJ. 



D. Negligible Coulomb energy 

Now, going back to the case of two qubits detected 
by two QPCs, a limit which is easy to analyze is when 
the rates do not depend on whether the doubly occupied 
state is involved or not, i.e. F^^^^ = F^^^. This regime 
is e. g. obtained if the Coulomb energy can be neglected 
U = Q. The average current is then 



pz pz 



(38) 



and the noise is the same as for a non-interacting quan- 
tum do1j2^ 



inr + (rij)^ 



The occupation probabilities of the island are 



(39) 



pa _ 
zz 



in? 



zz 



2riF|, 
(ri + rfj) 



(40) 



(Fi-f-Ff,)^ 

Now, looking at the symmetric situation where the left 
and right tunnelling rates are equal (F^ — Tr — F) 
and the qubits are equally strongly coupled to the SET 
{5Tl = STl = (5F), we have the following expressions for 
the time needed to distinguish the different qubit product 
states, 



AD 

ms 


T 


,BC 

''ms 


= CO 


AB 


,AC 

~ ''ms ~ 


,BD 

''ms ~ 


,CD 

ms 


ms 



6T' 



(41) 



to lowest order in ST/T. The states \A) = \ ||) and 
1^) = I TT) have the largest difference in average current, 
and are the fastest to distinguish. The states \B) = \ 
) and |C) = I ti) have the same average current, and 
thus we cannot tell them apart in this measurement. To 
distinguish all product states one has to measure for at 
least as long as the longest measurement time t^. The 
corresponding quantum efhciencies are 



AD 



„AB 



1, = 



0, 



CD 



(42) 



again to leading order in 5T /T. Here we may note that 
the quantum efficiency is unity only for telling the states 
\A) and |D) apart. The states \B) and |C) we can not tell 
apart, but they are still dephased, giving a zero quantum 
efficiency. 

For the other pairs of states, we find a quantum effi- 
ciency of one half, indicating that the qubits are getting 
entangled with a degree of freedom which is not mea- 
sured. This degree of freedom is the island charge state, 
and how the quantum efficiency can be improved by also 
detecting this will be discussed in Sec. |Vl 

In Refs. [26II27II28I , measurements which are unable to 
discern certain product states are used to create entan- 
glement. Our measurement setup can not be used for 
entangling the qubits, since it always projects the qubits 
on a product state. As will be shown in Sec. |Vl the in- 
ability to separate the \B) and |C) states is not inherent 
to the coupling between the qubits and the measurement 
device, but rather a consequence of measuring only the 
current. 



E. Large Coulomb energy 

A maybe more realistic situation is when the Coulomb 
charging energy U of the island is large, and especially 
larger than the applied drain-source voltage. The in- 
rates T I to the doubly occupied state are then zero for 
all qubit-states \z), and the doubly occupied state of the 
SET island will be unoccupied. In this situation the 
steady state island probabilities are 



pa _ 



2ri + F|j' 



Pt = — , P% = 0, 

2F|_+F|j' 



resulting in the average steady-state current 



pz pz 



and Fano factor— 



2Fi + F|,' 



/ 



4(ri)^ + (Ff,) 

(2Fi+F|j)2 



(43) 



(44) 



(45) 



We may note that the expressions for the current and 
Fano factor can formally be obtained from Eqs. ([55)) and 
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(139)1 by letting Tfj r|j/2. Again looking at the sym- 
metric situation where the left and right tunnelling rates 
are equal (F^ = Tji = F), and the qubits are equally 
strongly coupled to the SET {5Tl — STl = ST), we have 
the following expressions for the time needed to distin- 
guish the different qubit product states, 



t. 



AC 



,BC 



,BD 



5 F 
F 



,AB _ .CD 

ms ms 



15 F 



15 



ST' 



(46) 



where we again assumed weak coupling to the qubits 
^F <C F, and presented the results to leading order in 
ST/T. The states \A) = | H) and \D) = \ TT) again 
have the largest difference in average current, and are 
the fastest to distinguish. The states \B) = \ |t) and 
|C) — I ti) no longer give the same average current, and 
thus we can tell them apart. To distinguish all prod- 
uct states one has to measure for at least as long as the 
longest measurement time . The corresponding quan- 
tum efficiencies are 



„AD 
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^AC ^ ^BD 



10' 

1 

5' 



AB 



CD 



(47) 



to leading order in ST/T. Here all quantum efficiencies 
are finite, since we can distinguish all product states. 
Also, the quantum efficiency for separating the |^) and 
\D) state is lowered to 9/10 compared to unity in the 
noninteracting case, which is due to the differences in 
island state probabilities between the \A) and \D) state. 
In fact, by choosing T l = T ,T r = 2T ,ST l = ST and 
STfi = 2ST the quantum efficiency rj"^^ is again unity, 
and the measurement times and quantum efficiencies are 
given by the expressions for the noninteracting case in 
Eqs. gl]) and JMl)- 



both to a and c, so this waiting time is exponentially 
distributed with parameter F^ -I- Tr. By labeling the h 
state waiting times also with the destination state (a or 
c), one gets two separate waiting times (t^ and tJlf), each 
exponentially distributed, with parameter T r and F^ re- 
spectively. The average waiting time equals the inverse 
rate, e. g. (t° ) = F£^, and the variance is decreasing as 
F^^n^^, where n is the number of tunnel events (a — > b) 
where has been sampled. The average number of tun- 
nel events increases linearly with time, = iP°2Fi, 
Uba = tP^TR, Ubc = tP^T'^ and n, = tP^2F^. 

Now, including the coupling to the qubits, each of the 

four rates may have one of two possible values F^|^ or 

F^l^. The combined probability distribution for the four 
different waiting times will now depend on the qubits' 
state. We again use Eq. as the definition of the 
measurement time, with the difference that the proba- 
bility distribution now is a function of four continuous 
random variables. 

In the weak coupling regime (5F ^ F the individual 
distributions approach a Gaussian form, in the relevant 
regime of all Ua, n},a,nhc, nc ^ 1. The overlap for the 
distributions is 



dti^ \l Pz\ (t-w ; ^)Pz2 {^W ; ^) — ^ ^ ^ ^ 



(48) 



where is the average number of a ^ & tunnel events, 
neglecting the coupling to the qubits. This expression 
is valid to lowest non- vanishing order in all ST/T. Since 
the waiting times are independently distributed random 
variables we get the total overlap as a product of the four 
individual overlaps, giving the measurement time 



ms \ 



'S)l{zi,Z2) 



+ Qr{zi,Z2) 



Pr 



PbSTl 
2 Tr 



Pb jST'^ 
2 F'^ 



21 



Pr 



T' 



(49) 



V. MEASURING THE ISLAND STATE 

When the tunneling rates are slow enough, one may 
perform a time-resolved measurement of the island 
charge state, e.g. using another RF-SET^". Plotting 
the island charge as a function of time, each tunnel event 
appears as an almost vertical line, and from this informa- 
tion it is straightforward to determine the full counting 
statistics of the electrons passing the detector, as shown 
by Gustavsson et al^ 



A. Measuring the full time-trace 

The time the system waits in island state a (t,° ) and c 
(t^) are random variables, exponentially distributed with 
parameters F^ and F'^, respectively (neglecting the cou- 
pling to the qubits). ^From state 6, the system escapes 



where 6l/_r(2i,Z2) = 1 if the states \zi) and \z2) differ 
for the left/right qubit, and Ol/_r(-Z1j ^2) = otherwise. 
Comparing with the weak coupling expression of the total 
dephasing rate in Eq. (|29p , we find the quantum efficiency 
of recording all tunnel events to be unity, to lowest order 
in i5F/F. Thus, during the measurement, all the informa- 
tion about the qubits' state which is extracted is trans- 
ferred into information about the four different tunneling 
rates. 



VI. QUANTUM EFFICIENCY OF CURRENT 
MEASUREMENTS REVISITED 



Knowing that the full information is in the distribution 
of the rates, we can now reexamine the current measure- 
ment in terms of how much information it gives about 
the rates. In general the current is a function of all four 
rates, and the analysis gets complicated. For simplicity 
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we here only reanalyze the two cases we analyzed ana- 
lytically above: negligible and large Coulomb energy. In 
both these cases there is only one independent non-zero 
rate per junction. Thus the current is a function of only 
two independent random variables, and we may straight- 
forwardly visualize the situation. In Fig. [5] the pair of 
rates corresponding to the four different qubit product 
states A-D are shown as dots in a two-dimensional rate 
diagram. In the weak coupling regime the distributions 




' \ \ 

FIG. 5: A schematic representation of the different two- 
dimensional rate distributions corresponding to the 4 different 
qubit product states A-D. A current measurement is equiva- 
lent to projecting the distributions on the axis la in the case 
of negligible Coulomb energy, and on in the case of large 
Coulomb energy. 

are the bivariate Gaussians, centered around the corre- 
sponding state. The width of the distributions decrease 
with time as The green circles indicates where the 

distributions decreased to 1/e, at the measurement time 
^mf = = = ^mf • Siucc the distance between A 
and D is longer, these states separate earlier, and as 
is shown by Eq. (|49p this measurement time is only half, 

J.AD _ J.AB /o 
ms ms I 

The one-dimensional distribution of the current 
is given by integrating the two-dimensional rate- 
distribution along the axis perpendicular to the vector^ 

[9r,/(rL,Ffl),9r„/(rL,Ffl)]. (50) 

In the case of negligible Coulomb energy, the current 
formula is symmetric in the two rates, and the projec- 
tion axis is diagonal, as shown by the /a-axis. It is 
obvious that the separation of the states A and D has 
full quantum efficiency, since the line connecting the two 
states fall on the current axis. In general the quantum 



efficiency is given by rf^^'^ = cos^a^^^^, where a^^^^ is 
the angle between the current axis and the line connect- 
ing the states \z\) and |z2). Since the pairs of states 
AB, AC, BD, and CD are connected by horizontal or 
vertical lines, the corresponding quantum efficiency is 
cos^ (7r/4) = 1/2. Furthermore, the line connecting the 
B and C states is orthogonal to the current axis, giv- 
ing a zero (cos^ (''■/2) = 0) quantum efficiency, i.e. the 
projection of the two distributions always overlap com- 
pletely. This simple geometrical interpretation is in full 
agreement with Eq. ([1^ . 

In the other analytically tractable case of large 
Coulomb interaction, the current axis /f, forms an an- 
gle 7 = arctan(l/2) with the horizontal F^-axis. Noting 
that cos^7 = 4/5, and that cos^(7 — 7r/4) = 9/10 we 
see that the geometrical interpretation also reproduces 
Eq. g7|. 



A. Measuring average island charge 

Measuring the full time-trace of all tunnel events has 
so far only been demonstrated for tunnel rates on the or- 
der of kHz'^*'. For faster readout schemes a more realistic 
setup would be to measure the average island charge (n) , 
in addition to the average current through the device. 
The time-dependent distribution of the average charge 
(n) can also be analyzed in terms of the rate distribu- 
tions. One may note that the vector corresponding to 
the average charge measurement is perpendicular to the 
current measurement vector in the symmetric setup, i.e. 
F^ = Vji for negligible Coulomb energy, and F^ — 2Fi 
for large Coulomb energy. In this case the distribution 
of average charge is independent of the current distribu- 
tion, and we recover full quantum efficiency in separating 
all product states, even if the full time-trace cannot be 
recorded. 



VII. CONCLUSION AND DISCUSSION 

We have described how to read out two charge qubits 
using a single SET. After deriving general kinetic equa- 
tions for the whole system, we focussed on the properties 
of single-shot readout. Any interdot coupling during the 
measurement will reduce the fidelity, so we considered the 
case when the inter-dot qubit tunneling is turned off dur- 
ing measurement (A^, = = 0). Detecting the charge 
of the SET island in real-time corresponds to measuring 
the individual tunnel rates, and is always quantum effi- 
cient. The current is a function of the tunneling rates, 
and the quantum efficiency in this case depends on which 
qubit states should be separated, as well as the biasing 
conditions of the SET. In some symmetric cases, a com- 
bined measurement of average current and average island 
charge gives full quantum efficiency for separating all four 
product states. The setup can not be used for entangling 
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the qubits, since the measurement always projects onto 
a product state. 

The possibihty to have one readout device for two 
qubits can be used to optimize the design of many-qubit 
circuits. One drawback is that the qubits always have to 
be read out simultaneously. This may not be a severe 
problem, since in many quantum algorithms the readout 
is performed at the end. Also, there is an obvious advan- 
tage in reducing the number of readout Uncs attached to 
your circuit. This has to be weighed against the reduced 



signal-to-noise ratio in the readout, coming from the fact 
that four different states should be separated, instead of 
two. 
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